Line Widths of Single-Electron Tunneling Oscillations: 
Experiment and Numerical Simulations 

Jonas By lander, Tim Duty and Per Delsing 

Dept. of Microtechnology and Nanoscience, Chalmers University of Technology, SE-412 96 Goteborg, Sweden 

Abstract. We present experimental and numerical results from a real-time detection of time-correlated single-electron 
tunneling oscillations in a one-dimensional series array of small tunnel junctions. The electrons tunnel with a frequency 
f — I/e, where / is the current and e is the electron charge. Experimentally, we have connected a single-electron transistor to 
the last array island, and in this way measured currents from 5 fA to 1 pA by counting the single electrons. We find that the 
line width of the oscillation is proportional to the frequency /. The experimental data agrees well with numerical simulations. 



INTRODUCTION 

We have recently reported a measurement of a very small 
electrical current, /, (5 fA-1 pA) by direct counting of the 
single electrons that tunnel through a one-dimensional 
series array of metallic islands separated by small tun- 
nel junctions 1 1 ]. At low temperature, and when the array 
junction resistance is greater than the Klitzing resistance, 
/?k ~ 26 ki2, charge is localized and an excess electron 
charge, e, on one island polarizes the neighboring is- 
lands, thus repelling other electrons |2j,|3[. This poten- 
tial profile, extending over a distance of M — ^JCpjCo 
islands, is often called a "charge soli ton"; here C\ and 
Co are the junction capacitance and island stray capaci- 
tance, respectively. Electrostatic repulsion creates a lat- 
tice of solitons in the array, which above a certain thresh- 
old voltage starts to move as a whole, leading to an oscil- 
lation of the potential of any island inside the array with 
the mean frequency 



f = I/e 



(1) 



In (jj, the current was injected through a tunnel junc- 
tion from the array into the island of a single-electron 
transistor (SET) [|4j,|5|], a sensitive electrometer. The SET 
was embedded in a resonant LC circuit and operated in 
the radio-frequency mode (RF-SET) |6]. By monitoring 
the output signal from the RF-SET, we were able to de- 
tect the single-electron tunneling events in real time, and 
the power spectrum showed a peak at the frequency Q, 
thus demonstrating time correlation. 

This type of electron counter can potentially be used 
to measure small currents with very good accuracy, and 
without any need for calibration since the only parameter 
involved in is a natural constant. 



In this paper we compare the measured spectral line 
widths with those obtained from a computer simulation. 



EXPERIMENTAL TECHNIQUES 

We made the device out of Al/AlOx with an SET resis- 
tance of R SET — 30 ki2 and charging energy Eg /kg = 
1.6 K, where kg is Boltzmann's constant. We obtained 
a small-signal charge sensitivity of <5g = 2-10~ 5 e/Hz 1 / 2 
and a bandwidth of 10 MHz at an RF carrier frequency of 
358 MHz. The N = 50 junction array, with junction pa- 
rameters =940 kQ., C A =0.42 fF and C =0.030 fF, 

had the charging energy £^ lay /kg = 2.2 K per junction, 
and a soliton size of M ks 3.7 islands (that is, N 3> M). 

We performed the experiments in a dilution refrigera- 
tor at 0.03 K and a parallel magnetic field of 475 mT. The 
device was superconducting, but with a suppressed gap, 
A/ kg « 0.6 K. At this field, the Josephson coupling en- 
ergy is very small in our junctions: Ej /kg < 10 mK <C T. 
Therefore Cooper pair tunneling is strongly suppressed, 
and we see only quasiparticles. Furthermore, we take ad- 
vantage of the subgap resistance, R§q 
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for an electron temperature T e « 0.15 K, es- 



timated using parameters for electron-phonon coupling 
and Kapitza resistance 17!]. This makes the onset of cur- 
rent more gradual than in the normal state, thus reducing 
the sensitivity to bias voltage fluctuations. 



NUMERICAL METHODS 



We have numerically simulated the electron transport in 
the array using a direct Monte Carlo method jj, |3[ . We 




FIGURE 1. Experimental (a) and simulated (b) power spectral densities (solid lines) with Lorentzian fits (dotted). The currents 
were 79, 220 and 370 fA (a), and 71, 221, 374 fA (b). Insets: Time series for / = 80 fA (in (b) before filtering), where each peak 
represents one electron passing through the SET. Comparison of fitted frequencies, /o, (c) and relative line widths, y//o, (d) for 
experimental (dots) and simulated (circles) spectra. Solid line in (c): / = ef. Inset in (c-d): Array /— V curve at By, = 475 mT. 



calculate the tunnel rates in each junction by using the 
"orthodox" theory of single electron tunneling 1 8], where 
the probability of a tunneling event per unit time is fully 
determined by the change in free energy of the system. 
The time step between such events is determined by a 
random variable that mimics the stochastic nature of tun- 
neling. Below the superconducting gap, V < 2A/e, we 
have phenomenologically introduced the subgap resis- 
tance, R§g- For simplicity, however, we disregard the su- 
perconducting density of states above the gap. The simu- 
lations are idealized in the way that we disregard disorder 
and external noise, as well as bias fluctuations or drift. 

We model the charge detection by sampling the poten- 
tial of the last array island at a rate corresponding to the 
bandwidth of the LC circuit. 



few fA during the measurement time. The deviation from 
Q in Fig-HJc) for / > 400 fA arises probably because of 
an increasing asymmetry of the peak, which causes the 
fit to underestimate the frequency. 

In conclusion, we have shown that the line width of 
the single-electron tunneling oscillations is proportional 
to the frequency, and that the line width extracted from 
our simulations agrees well with the experimental data. 
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RESULTS AND DISCUSSION 

From the experimental data we have subtracted a back- 
ground to account for amplifier and external noise, in- 
cluding an approximate 1 // term. Then, both experimen- 
tal and numerical power spectra are fitted to the function, 

S(f)=af/((f-f ) 2 + f), (2) 

a Lorentzian around a center frequency /o, with a half 
width 7, see Fig.^a-b). We find that /o agrees with Q, 
see Fig. Qc), and that y is proportional to the current: 
Y/fo w 0.4, see Fig.Ed). 

For / < 200 fA, the measured peak is broad compared 
to the simulation, likely due to the difficulty in maintain- 
ing a stable bias for very small currents. We estimate the 
noise from the biasing circuitry and thermal EMF to be 
of the order of 10 jJ.W, causing a current fluctuation of a 
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